The three site interacting spin chain in staggered field: Fidelity vs Loschmidt echo 
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We study the fidelity and Loschmidt echo of a three site interacting XX chain in presence of 
staggered field which consists of special types of quantum phase transition points due to change in 
the topology of the Fermi surface, apart from normal phase transition points from gapped to gapless 
phases. We focus on the effect of such types of phase transition points on the fidelity, which is the 
static part of the problem, and Loschmidt echo of this spin chain. We find that on one hand, fidelity 
is able to detect only the boundaries of the gapped to gapless phase transition and not the phase 
transition from two Fermi points to four Fermi points which lie within this gapless region. On the 
other hand, Loschmidt echo is just able to detect one special point in the entire phase diagram. We 
present the arguments for such a behavior. 

PACS numbers: 



I. INTRODUCTION 

The desire of building a quantum computer to solve 
quantum problems efficiently has lead to the development 
of quantum information theory and other related fields. 
Fidelity, which is basically an overlap of the ground state 
wave function at two slightly different values of the pa- 
rameter, is an important component of the quantum in- 
formation theoryir^. It is important because it can de- 
tect a quantum critical point without knowing anything 
about the order parameter of the system, which other- 
wise is the conventional way of looking at the quantum 
phase transition (QPT)^£. Fidelity has been successful 
in detecting various types of quantum phase transition 
points, for example, ordinary critical points separating 
two gapped phases through a gapless poinfc^ or topo- 
logical QPTs£ or QPT in Bosc Hubbard modeli, but its 
usefulness in a general scenario is not very clear. We 
present a model where the fidelity detects the onset of 
quantum critical region but not the two different phases 
within this quantum critical region. 

One of the other important quantities in quantum in- 
formation theory is Loschmidt echo (LE)£~— . LE is the 
measure of the overlap of the same initial state, the 
ground state of the initial Hamiltonian H(X), evolving 
under the influence of the two Hamiltonians, H (A) and 
H(X + 5). This overlap shows a dip at the quantum crit- 
ical point, thus enabling us to detect a quantum critical 
point. In the language of quantum information theory, 
LE can be used to measure the quantum to classical tran- 
sition of a qubit coupled to an environment making it 
an important component of quantum information theory. 
LE was actually introduced in connection to the quan- 
tum to classical transition in quantum chaosii - — and 
now extended to various other systems like Ising model^, 
Bose-Einstein condensate models and Dicke modet^. It 
has also been studied experimentally using NMRji£~— . 
Once again, LE has been used to detect the quantum 
critical points but its generality is yet to be explored. 

In this paper, we point out the inability of fidelity or 
LE to detect certain special class of critical points by tak- 
ing the example of staggered transverse field in a three 



spin interacting spin-1/2 XX chain. This study is an im- 
portant addition to the literature as no one has yet looked 
at the fidelity and LE of a Hamiltonian which also un- 
dergoes a very unique type of quantum phase transition, 
namely, from two Fermi points to four Fermi points. It is 
not apriori obvious the effect of such points on the fidelity 
and LE of the model. It is to be noted that many interest- 
ing experimental observations, like magnetic properties 
of solid 3 He, were explained by the inclusion of multispin 
interaction in the Hamiltonian^i. 

The outline of the paper is as follows: The model is 
described briefly in Section II. Section III discusses the 
statics, i.e., the ground state fidelity of the model and 
the dynamics of ground state Loschmidt echo is studied 
in Section IV. We present our conclusions in Section V. 



II. MODEL 

In this section, we quickly go through the ground state 
phase diagram of the three spin interaction Hamilto- 
nian in presence of a staggered field h s given by the 
Hamiltonian^ 
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Performing standard Jordan Wigner 

Fcrmionizatio n 23 i 24 , we get 



H ts = -77 X^4 c J+ 1+ S t +i c j) + T^( c I c i+ 2 ~ r 



C-A_\nC"i 



(-l) j h s (c] Cj \) (2) 

To diagonalize the above Hamiltonian, it is useful to in- 
troduce two types of spinless fermions on the odd and 
even sub-lattices as shown below: 

C2n-i = a„-i/2 an d c 2n = b n . 



Substituting this in Eq. [2] and Fourier transforming it 
gives, 



Hts — /_^ Hk 



k 



where 



e a (k) 

e b (k) 
e a b(k) 



e a (k)a\a k + e b {k)b\b k + e ab (k)(alb k + b\ai$\ 



= — cos(fc) + I 

= y cos 0) - ' 
= -Jcos(fc/2) 



(4) 

It is easier to understand the fidelity and LE calcula- 
tions if we rewrite the Hamiltonian H k in the four basis 
described by 1 0,0), \a k ,0), |0,6fc), \a k , b k ) where the first 
index represents the presence or absence of a k — particle 
and the second index corresponds to that of b k — particle. 
In these basis, H k is given by 



H k = 
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e b (k) 



(5) 



It is clear from the above representation that to di- 
agonalize the above Hamiltonian, only the basis \a k ,0) 
and |0, b k ) needs to be rotated since the Hamiltonian in 
the other two basis is already diagonal, i.e., there is no 
mixing along these two directions. Let us denote the two 
new directions for the two new quasi-particlcs a k and f3 k 
as |afe,0) and \0,(3 k ) which diagonalize the total Hamil- 
tonian giving four eigen energies 0, E—(k), E + (k) and 
E_(k)+E + (k) corresponding to the four directions |0, 0), 
|afc,0), |0,/3fc) and \a k ,/3 k ), respectively. These two new 
directions are given by 



|afc,0) = cos — \a k ,0) -sin — 10, b k ) 
|0,/3 fe ) = sin— \a k ,0) + cos — \0,b k ) 



where 



t&ii9 k (h s 



— cosfc/2 



(0) 



(7) 



with J set equal to unity through out the paper, and 
their energies are 



E± (k) = ^ cos k ± y/h 2 s + cos 2 (fc/2) . 



(8) 



Hence, the Hamiltonian is now given by 

Hts = Y, E_{k)a\a k + E + {k)$[p k . 



E 

k 



The ground state of this Hamiltonian corresponds to all 
the modes with negative energies filled and hence depend- 
ing upon the values of the Hamiltonian parameter, the 
system has various phases. We discuss briefly the various 
phases below: 




FIG. 1: The phase diagram of the three spin model in presence 
of staggered field. The solid lines correspond to h s = ±Js,/2 

critical line and the dotted lines correspond to h s = ±\ 



• When h s > J 3 /2: E + (k) > and £L(fc) < 
for all k, and hence the ground state for each 
mode is |afc,0) with total ground state energy 
Eg = Y^ k E-{k). This is the Antiferromagnetic 
(AF) phase. 



When spffl - 1 < \h\ < J 3 /2: Some of the k 
modes from the E + (k) branch become negative 
whereas E-(k) is negative for all the modes. This 
phase has two Fermi points arising due to zeros of 
the E + (k) branch. In this region, the ground state 
for a given mode can be |a/c,0) or \a k ,/3 k ) depend- 
ing upon whether E + {k) branch is empty or filled, 
respectively. The total ground state energy is given 

by 
J2E-(k)@(-E_(k))+E + (k)e(-E + (k)), 



where 9 is the Heaviside function. We call this 
phase as spin liquid I (SLI) phase. 



• When < \h\ < ^/Jf/4- 1: In this limit, E_(k) 
also crosses zero for some modes and hence there 
are three possible ground states for a given mode 
k depending upon the signs of the energies E± (fc) 
given by \a k , 0), |0, 0) and \a k , (i k ). Once again the 
ground state energy is the sum over all the modes 
with negative energies for each branch as written 
above. This phase is called spin liquid II (SLII) 
phase. 

Following the above arguments, the phase diagram of the 
model is shown in Fig. [TJ We now try to capture these 
phase transitions using the fidelity approach, especially 
the phase transition from SLI to SLII and later discuss 
them in the light of Loschmidt echo. 



III. FIDELITY 

As mentioned in the introduction, the ground state 
fidelity is defined as the overlap between the two ground 
state wave functions at slightly different parameter values 
and is given by 

F = (y G (h s )\V G (h 3 +S)) = Y[($ h (h,)\$ k (h.+S)) =Y[F k 

k k 

(9) 

where \^c(hs)) is the total ground state for the given 
parameter value and |$fc(/i s )) is the ground state for the 
k — th mode. Thus, while evaluating F k , one has to care- 
fully identify the ground state for the k — th mode de- 
pending upon the sign of E±(k) . The various possibilities 
are: 

. F k = (a k ,0\a k ,0) s = C o S ( MM^±g2 ) 

• F k = (a k ,p k \a k ,p k ) s = (0,0|0,0) 5 = 1 

• F k = (afc,0|afe,/3fe}(5 = (afc,0|0, 0)«5 = with simi- 
lar cross products also equal to zero. 

Here, we use the notation of the bra/ket without any 
subscript for the field h s and with a subscript S for the 
field h 8 +5, which will be followed through out this paper. 
We have also used the orthogonality of the basis states in 
the above steps. It is to be noted that within the entire 
gaplcss region including SLI and SLII phase, there is at 
least one mode for which F k — and hence fidelity in the 
entire gapless region is zero. This is also shown in Fig. 
® 
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FIG. 2: Variation of fidelity as a function of the staggered 
field h s for two different values of J3. The SLII phase exists 
only for J3 = 4 but this phase is not captured by the fidelity. 
Fidelity only detects the gapped to gapless phase transition 
of the model which occurs at h s = J3/2. 



IV. LOSCHMIDT ECHO 

The ground state LE, defined as the overlap of the ini- 
tial wavcfunction, which is the ground state of the initial 



Hamiltonian, evolving under two slightly different Hamil- 
tonians, is given as 



LE=\(*(h s ,t)\y(h s + 6,t))f 



(10) 



where |#(ft 8 ,i)) = e iff ^)*|* G ) and \^(h s + S,t)) = 
e iH(h s +8)t^ G ^ |^, G ^ bdrig the ground state for the 

Hamiltonian H = H(h s ). In the momentum represen- 
tation, as in the case of fidelity, one has to write an ex- 
pression for the Loschmidt echo by taking the correct 
ground state for each mode. Thus, 

LE = H\(^ k (h s )\e iH ^' +s ^k(h s )}\ 2 

k 

= \{LE k (11) 

fe 

Let us calculate LE^ for different possible ground states. 

• if |$ fc (/i s )) = K,o) 

2 



LE, 



k 



(a k ,0\e iH ^+ 5 ^\a k ,0) 



(12) 



Since |afc,0) is the eigenstate of H(h s ) and not 
H(h s +S), we need to rewrite it in terms of \a k , 0)s- 
Using Eq. [51 we find that 



\a k ,0) = cosr) k \a k ,0) s - sinri k \0,/3 k )s, 



(13) 



where 2rj k = 9 k {h s ) — 9 k (h s + S). Substituting the 
above transformation in Eq. 1121 we get 



LE k = 1 - sin 2 (2r/ fc ) sin 2 f AE ^ 



(14) 



where AE(k) = E+(h s + 6,k) - E_(h s + 5,k) = 
2sJ(h s + S) 2 + cos 2 (fc/2). Clearly, the three spin 
term has nothing to do with in AE and does not 
influence LE^ though it will affect the magnitude 
of the dip in LE as shown in Fig. 3. This is be- 
cause the number of k— modes which influence the 
Loschmidt echo changes as J3 changes. 

• If \& k (h.)) = |0,0) or \a k ,p k ), then LE k = 1 as 
these are also the eigen vectors of the Hamiltonian 
H k (h s + S). 

From the above two points, it is clear that the LE will 
show a dip when AE is minimum at which sin 2r\ k is non- 
zero. On the other hand, the ground state at h s + 6 = 
for the mode k = 7r at which AE = is \a k ,(3 k ) for 
J3 7^ and hence LEfc =7r is unity making this critical 
mode irrelevant in the dynamics. It is the other modes 
for which E + (k) > and £L (fc) < which actually influ- 
ence the behavior of the LE and is one of the interesting 
observations of the paper. Since the critical mode is not 
involved in the dynamics, various power-law scalings of 
LE are not present in this model. There is another way 
of looking at this point: Consider Eq. [5] As mentioned 
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FIG. 3: Variation of the LE with the staggered field h 3 
different values of three spin interaction term. 



for 



before, |0, 0) and \a,k, bk) do not mix and we can concen- 
trate on |afc,0) and |0, bk) basis. In these two basis, H^ 
can be written as 



j 3 

— cos kl + h s a z — cos(k/2)a 0L 



(15) 



where second and third term combined together resem- 
bles the Landau -Zener Hamiltonian, though without the 
linear time dependence 2 ^. Clearly, the dynamics is gov- 
erned by the minimum energy gap corresponding to the 
Landau Zener Hamiltonian 2 ^ which occurs at h s = 
and hence only this point is detected by the LE. We call 
h s = and k^ = ir as the dynamical critical point and 
dynamical critical mode, respectively. Such an observa- 
tion was also reported in the context of defect genera- 
tion for a similar model where the authors claimed that 
the identity term only changes the phase of the evolv- 
ing wave function and hence does not contribute to the 
defect generation 2 ^ 

We also studied the dependence of LE on system size 
and the coupling strength d at the dynamical critical 
point h s + 6 = 0. The LE decays exponentially with 
the system size N whereas it decays exponentially with 
S 2 for both the regions, \h s \ > J3/2 and \h s \ < J3/2 as 
shown in Fig. [4] To understand this behavior, let us 
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FIG. 4: Variation of log LE with system size N and the cou- 
pling strength S for h s + S — 0. 



concentrate on S =log LE which is given by 

S = £log(l-sin 2 2 % sm 2 ^- r 
k ^ 



(16) 



As mentioned before, AEk never vanishes as the critical 
mode does not contribute to the decay of the LE. Since 
sin 2 2rjk is very small as it is the difference between two 
approximately equal angles (critical mode is not involved 
in the calculation), we can write 



S 



£4,Mn 2 (^) 



E 



(5 2 sin 2 (tcosfc/2) 
cos 2 (fc/2) 



(17) 



which explains the 5 2 dependence of S. It is difficult to 
understand the system size dependence of the LE since 
the critical mode does not contribute to the decay and 
hence no further simplifications can be done. But we 
do know that the mode closest to the dynamical critical 
mode fcj? = 7r for which E-(k) < and E + (k) > con- 
tribute maximum. Expanding Eq. [17] around k = tt for 
small time i, wc find that 



E 



S 2 t 2 oc N 



and is shown in Fig. 2J 

Though we considered a spin-1/2 model, such type 
of excitation spectrum, i.e., identity matrix and a Lan- 
dau Zener type Hamiltonian, also exists in hard core 
Bosc-Hubbard model in the presence of a period two 
supcrlattice21. We expect similar results to hold even 
in this hard core boson model which has various phases 
like superfluid, Mott insulator, hole vacuum and particle 
vacuum phase. 

At the end, we note that there are a couple of works 
studying the fidelity and LE of models containing multi- 
spin interaction in Ising type Hamiltonian a 28 ' 29 . In Ref. 
|28| . the effect of some complicated three spin interaction 
on LE is studied where the authors reported that a partic- 
ular term does not effect the position of the dip in LE and 
changes only the sharpness of the decay. On the other 
hand, we have studied a completely different model with 
special types of phase transition points and presented ar- 
guments on why LE is not able to detect the gapped to 
gaplcss phase transition points. Also, the physics of our 
model is very different from the model studied in Rcf. 
[28l due to the presence of staggered field. As a result of 
this, our model contains two types of quasiparticlcs mak- 
ing the handling of the problem very different from the 
other models studied till now in the context of fidelity or 
Loschmidt echo. 



V. CONCLUSIONS 

Usually, there is a sharp dip in the ground state fidelity 
or equivalently in the LE close to a quantum critical 
point. We have shown here an example of a unique quan- 
tum phase transition where neither fidelity nor LE can 
capture the entire phase diagram. While fidelity shows a 



dip at the boundary between the gapped to gaplcss phase 
transition and is unable to detect the two Fermi point to 
four Fermi point phase transition within the gapless re- 
gion, the Loschmidt echo shows a completely different 
picture. It is able to detect only one special point in the 
entire phase diagram which is the point corresponding to 
the minimum of the Landau Zener type Hamiltonian and 
can not capture the critical points of the total Hamilto- 



nian. This is because the dynamics is entirely governed 
by the Landau-Zcner type Hamiltonian and the identity 
matrix containing the three spin term only adds a phase 
to the wave function evolution. 
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